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Abstract 

PAMELA and ATIC recently reported excesses in cosmic rays. 
Since the interpretation in terms of DM annihilations was found to be 
not easily compatible with constraints from photon observations, we 
consider the DM decay hypothesis and find that it can explain the 
excesses compatibly with all constraints, and can be tested by dedi- 
cated HESS observations of the Galactic Ridge. ATIC data indicate 
a DM mass of about 2 TeV: this mass naturally implies the observed 
DM abundance relative to ordinary matter if DM is a quasi-stable 
composite particle with a baryon-like matter asymmetry. Technicolor 
naturally yields these type of candidates. 



1 Introduction 



The recent observations of Cosmic Ray spectra by PAMELA [1] and ATIC [2] indicate an 
excess, compatibly with other recent results [3, 4]. Unless it is due to pulsars or some other 
astrophysical source, this excess could be the first non-gravitational manifestation of Dark 
Matter. However, the interpretation in terms of DM annihilations into SM particles is possible 
only for a restricted class of DM models [5, 6, 7, 8] leading to an unobserved excess of gamma 
or radio photons, unless the DM density profile is significantly less steep than the NFW profile 
at galactocentric radii below 100 pc [9]. 

In this paper we show that the interpretation in terms of DM decays [10] circumvents these 
issues. In section 2 we compute the 'halo functions' that encode the astrophysical information 
relevant for understanding the energy spectra of and p produced by DM decays. In section 3 
we propose a model-independent interpretation of the PAMELA and ATIC observations in 
terms of DM decays, taking into account constraints from p data. In section 4 we show that 
the models that can explain the PAMELA and ATIC excesses are compatible with photon 
observations at gamma and radio frequencies, and single out the future observations capable 
of testing the decaying DM interpretation. We will find that: 

• When interpreted in terms of decaying DM, PAMELA e"*" and ATIC e"*" -|- e~ data can 
be accounted for by a DM particle with mass around 2 TeV and lifetime ~ 10^^ s. Decay 
modes yielding hadrons are constrained by PAMELA p data and should not exceed 10%. 

In section 5 we explore the possible theoretical realizations of such scenario, finding that: 

• If DM particles carry an anomalous global charge B', the observed DM abundance, 
^dm/^b ~ 5, is naturally obtained for Mdm ~ 2 TeV. 

• The above situation can be realized assuming that DM is a composite state, suggesting a 
connection with technicolor: if the hypothetical DM is the lightest state carrying techni- 
baryon number B', a naive rescaling of QCD suggests a Mdm ~ 2 TeV mass. 

• If B' is violated by dimension-6 higher dimensional operators suppressed by some 'GUT 
scale' ~ 10^^ GeV, DM is quasi-stable and the desired DM life-time is naturally obtained. 

Section 6 summarizes our results. 

2 Astrophysics 

The computation of the fluxes of SM particles k = {e~^,p, . . .} from DM decays is analogous to 
the case of DM annihilations: in brief one needs to replace the "annihilations" source term 
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Figure 1: Left: The uncertain 'halo function' I {Xd) of eq. (5) that encodes the astrophysics of 
DM decays into positrons and their propagation to the Earth. The diffusion length is related 
to energy losses as in eq. (6). Right: The p astrophysical function R(T) of eq. (10), computed 
under different assumptions. In both cases, the dashed (solid) [dotted] bands assumes the min 
(med) [max] propagation configuration of eq. (4) and eq. (9) respectively. Each band contains 3 
lines, that correspond to the isothermal (red lower lines), NFW (blue middle lines) and Moore 
(green upper lines) DM density profiles. 



where dNk/dE is the spectrum of particles k produced by one decay or annihilation, p{x) the 
DM energy density at x and Mdm its mass. Furthermore DM sub-halos can enhance the DM 
annihilation signals by a potentially large 'boost factor' Bk > 1; for DM decays one expects a 
boost factor negligibly different from unity. 

We briefly summarize aspects of and p^ in our galaxy in order to present the astrophysical 
functions, plotted in fig. 1, that allow to compute their fluxes at Earth in the generic case. 

2.1 Positron propagation 

The flux per unit energy of ultra-relativistic positrons is given by $e+(^; x,E) = f /Air where the 
positron number density per unit energy, f(t, x, E) = dN^+/dE, obeys the stationary diffusion- 
loss equation: 

-K{E).V'f-^{b{E)f)=Q^ (3) 

with diffusion coefficient K{E) = Kq{E / GeYf and energy loss coefficient b{E) = E"^ / ( GeV • 
te) with te = 10^^ s. They respectively describe transport through the turbulent magnetic 
fields and energy loss due to synchrotron radiation and inverse Compton scattering on galactic 
photons. Eq. (3) is solved in a diffusive region with the shape of a solid flat cylinder that 
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sandwiches the galactic plane, with height 2L in the z direction and radius R = 20 kpc in the 
r direction. The location of the solar system corresponds to x = (r©, Zq) = ((8.5 ± 0.5) kpc, 0). 
The boundary conditions impose that the positron density / vanishes on the surface of the 
diffusive cylinder, outside of which turbulent magnetic fields can be neglected so that positrons 
freely propagate and escape. The values of the propagation parameters 6, Kq and L are deduced 
from a variety of cosmic ray data and modelizations. We adopt the sets discussed in [1 1] : 
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Analogously to the DM annihilation case [6, 11, 12], the solution for the positron flux at 
Earth can be written in a useful semi-analytical form 

^E'^-'iME^E-)) (5) 

where X£){E, E') is the diffusion length from energy E' to energy E: 

Xl = AKoTE 

the adimensional 'halo function for DM decays' I{Xd) fully encodes the galactic astrophysics 
and is independent on the particle physics model. Its possible shapes are plotted in fig. 1 
for the set of DM density profiles (isothermal, Navarro, Frank and White (NFW) [13] and 
Moore [14]) and propagation models that we consider. We see that, unlike in the case of 
DM annihilations, / is negligibly affected by the uncertainty in the DM density profile, and 
DM decays in the center of the galaxy (at 8.5 kpc from us) are never dominant. 
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2.2 Antiproton propagation 

The propagation of anti-protons through the galaxy is described by a diffusion equation anal- 
ogous to the one for positrons. Again, the number density of anti-protons per unit energy 
f{t, X, T) = dNp/dT vanishes on the surface of the cylinder at z = ±L and r = R. T = E — rup 
is the p kinetic energy, conveniently used instead of the total energy E. Since rrip ^ nie we can 
neglect the energy loss term, and the diffusion equation for / is 

- K{T) ■ V V + ^ (sign(z) / Konv) = Qp-2h S{z) T^^J (7) 

where: 

- The diffusion term can again be written as K(T) = Kq[3 {p/ GeV)*^, where p and [3 are 
the antiproton momentum and velocity. 5 and Kq are given in table (9). 

- The Vcomr term corresponds to a convective wind, assumed to be constant and directed 
outward from the galactic plane, that tends to push away p with energy T < 10 mp. Its 
value is given in (9). 
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- The last term in eq. (7) describes the annihilations of p on interstellar protons in the 
galactic plane (with a thickness oih = 0.1 kpc <^ L) with rate Fann = (^H+4^^^?T,He)c"p^'^Wp5 
where nn ~ 1/cm^ is the hydrogen density, nne ~ 0.07 nn is the Helium density (the factor 
4^/^ accounting for the different geometrical cross section in an effective way) and app^ is 
given by [12] 

_ / 661 (1 + 0.0115 r-°-^^4 - 0.984 mbarn, for T < 15.5 GeV , . 

^PP ~ I 36 mbarn, for T> 15.5 GeV ' 

- We neglect "tertiary p" , i.e. non-annihilating p interactions on the matter in the galactic 
disk. 

The set of propagation parameters that we adopt in the case for anti-protons has been 
deduced in [15]: 
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The solution for the antiproton flux at the position of the Earth $p(T, Tq) = fp/(47r)/ can 
be written as 

where pQ = p{rQ)- The 'halo function for DM decays' R{T) encodes all the astrophysics and is 
plotted in fig. 1 for the halo and propagation models that we consider. It depends negligibly on 
the DM density profile, unlike for the analogous function relevant for DM annihilations, plotted 
in fig. 4 of [6]. 



3 The PAMELA and ATIC data as DM decays 

With the goal of performing a model-independent analysis, we delineate the generic features 
implied by known physics. We recall that non-relativistic DM annihilations are equivalent to 
the decay of the s-wave DM DM two-body state with mass 2Mdm, that can only have spin 0, 
1 or 2, if DM is a weakly-interacting particle [5]. 

The decay of a DM is less constrained: since DM could be coupled only gravitationally to 
SM particles, its spin can be anything: 0, 1/2, 1, 3/2, 2, etc. Furthermore, since DM decay must 
be slow, decays into two SM particles do not need to dominate over multi-body decays; decay 
rates suppressed by helicity factors such as (me/MoM)^ are still phenomenologically interesting; 
DM decays might involve new particles (e.g. gravitino decays into gluinos and quarks). 

At the light of the results of [5] that single out decays into hard leptons as phenomenologi- 
cally promising we will study the following DM decay modes. For boson DM: 

DM e+e- , fi+ p- ,T+T- ,W+W- ,tt,hh (11) 
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Figure 2: A fit of the DM decays indicated in the legend to the PAMELA positron data. 



where the tt and hh modes are included as examples of hadronic modes. For fermion DM we 
consider: 

BM-^W^e^,W^fi^,W^T^. (12) 

3 body decays, such as DM i^i'u are also possible, but cannot be computed in a model- 
independent way. 

We fit the PAMELA and ATIC data taking into account as described in [5] the uncertainties 
on the astrophysical backgrounds, on the e^and p propagation, on the DM density profile, on 
the experimental data. 

3.1 PAMELA e+ data 

We start including only the PAMELA data about the positron fraction [1], that exhibit an 
excess above 10 GeV. Fig. 2 shows the quality of the best-fit as function of the DM mass for 
the different decay modes we consider, and fig. 3 shows the corresponding best-fit values of the 
DM life-time, that is r ~ 10^^ sec. The left column of fig. 5 show three possible fits. 

3.2 PAMELA e+ and e+ + e data 

Next, we explore the combined fit of PAMELA e"*" and the + e~ data, as measured by balloon 
experiments [2] and at higher energies by HESS [4]. We plot and fit the HESS data combining in 
quadrature a systematic ±20% systematic error with a ±15% energy-scale uncertainty with the 
statistical uncertainty on each data point [4]. Fig. 4 shows our results: the peak possibly present 
in the e"*" ±e~ energy spectrum, as measured by ATIC-2 and PPB-BETS, strongly restricts the 
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DM mass in GeV 



Figure 3: Best-fit values of the DM life-time suggested by the PAMELA excess, for the DM 
decay modes indicated on the legend. 



DM mass and the DM decay modes. Like in the DM annihilation case, the ijl^ and r^r^ 
modes provide good fits, now for a value of the DM mass twice as large: M TeV. Unlike 
in the DM annihilation case, the e"'"e~ mode does not provide a good fit, due to the following 
detailed argument, that can be weakened invoking additional astrophysical uncertainties. The 
problem is that the line remains too sharp. Indeed DM decays (unlike DM annihilations) 
are not strongly enhanced close to the galactic center, so that such e+ produced far from us 
and that therefore reach us at the price of significant energy losses are now a minor component. 
Furthermore, the new W^t^ decay modes possible for fermion DM also provide good fits for 
and e"*" + e~ data for any ^. A similar result holds for the various 3-body decays DM — t/, 
that we cannot compute in a model-independent way. 

3.3 PAMELA e+ and p data 

Fig. 6 shows the result of the fit of PAMELA and p PAMELA observations. The p data, that 
do not show an excess with respect to the expected astrophysical background, provide significant 
constraints, disfavoring decay modes that lead to hadrons. The situation with the DM — >■ 
W~^W~ ,ti, hh modes is similar to the DM annihilation case [5]: such modes are disfavored 
unless M > 10 TeV. We see that the same conclusion applies to the new semi-leptonic decay 
modes W^i^: choosing M ~ 2 TeV in order to fit the -\-e~ peak in ATIC-2 data is disfavored 
by PAMELA p data. This is illustrated in the upper row of fig. 5, where we summed over 
all leptons £ = {e, /i, r} as suggested by gauge invariance. We recall that we assumed equal 
propagation models for e"*" and p, as controlled by the thickness L of the cylinder where turbulent 



7 



85 




75 - 



'>t 70 



55 



300 



1000 



3000 



10000 



DM mass in GeV 



Figure 4: Combined fit of positron data from PAMELA and of + e data from ATIC, BBP- 
BETS, EC, HESS. 



magnetic fields lead to diffusion of charged particles. Our conclusions can be weakened assuming 
different propagation models for p and e"*". Unlike in the DM annihilation case it is not possible 
to weaken these conclusions by assuming a boost factor Bp <^ Bg. 

Summarizing, a good fit to all e^, e+ + e~ and p CR data presently available in terms 
of decaying DM requires decays into fi'^^" or r+r^, which are possible if DM is a boson. A 
fermion DM can decay into i'^i'u, which can also provide a good global fit, although we cannot 
compute the resulting energy spectra in a model-independent way. If decaying DM is only 
a fraction / of all DM, our results still hold up to the rescaling t ^ r/f. 

4 Photons: bounds and signals 

DM decays unavoidably lead to photons: at 7 frequencies due to brehmstralung of charged 
particles, and at radio frequencies due to synchrotron radiation emitted by in the galactic 
magnetic fields. Here we compare the flux of photons produced by DM decays with the existing 
bounds and we comment on future perspectives. There are three main robust constraints and 
we find that, unlike in the case of DM annihilations [9], they are satisfied even in the case of a 
NFW DM density profile. 
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DM with M = 1 TeV that decays into W-r 




Positron energy in GeV Energy in GeV p kinetic energy in GeV 




DM with M =A TeV that decays into t'^t 




Positron energy in GeV Energy in GeV p kinetic energy in GeV 



Figure 5: Examples of fits of (left), e+ + e~ (center), p (right) CR data, for DM decay 
modes into W^t^ (that leads to an unseen p excess), e"*"e^ (that leads to a too sharp peak), 
T^T~ (that leads to peak less sharp than suggested by ATIC data). A combination of the last 
two possibilities, and the intermediate DM — > decay gives the optimal fit. 
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Figure 6: Combined fit of different DM decay channels to the PAMELA positron and PAMELA 
anti-proton data, assuming equal propagation models for and p. 

4.1 HESS observation of the Galactic Center 

The differential ffiix of pliotons received from a given angular direction dfl is 



where dN^/dE is the photon spectrum produced in one DM decay: we here include only the 
model-independent contribution due to brehmstralung of charged particles. We conservatively 
impose that such contribution does not exceed at 3cr the observed flux in any data-point. The 
adimensional quantity J describes the uncertain astrophysics. Observing an angular region 
Ar2 = 10^^ centered on the GC, we find 

GC J for DM decays = {5.75, 28.9, 45.3} for the {isoT, NFW, Moore} profiles (14) 

for the quantity J defined by J ■ Afi = / dVl ■ dJ/dQ. In the case of DM annihilations, one 
instead has J ■ AQ = J dQ J^p/pQ^ds/rQ and, the analogous factor J acquires much larger 
values, such as J = 14700 for the NFW profile, implying that the PAMELA anomaly is not 
compatible with bounds from 7 observations, unless DM is distributed with an isothermal-like 
profile and/or the signal is enhanced by a large boost factor > 10^ ^ [9]. 

On the contrary this problematic issue is not present for the DM decay interpretation of 
the PAMELA excess, as shown by the blue continuous lines in fig.s 7. 
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4.2 HESS observation of the Galactic Ridge 



It is unclear which observation region gives the maximal sensitivity to 7 rays from DM decays. 
Observing the Galactic Center (GC) maximizes the signal rate, especially for cusped profiles 
such as NFW, where the DM density p grows as 1/r as r — 0, and especially for DM anni- 
hilations, where the signal rate is proportional to p^. However, the GC also has the highest 
background rate, as HESS observations suggest that it is polluted by astrophysical sources of 
7 rays. 

For DM decays, the signal is proportional to p rather than to p^, so the optimal strategy 
might be instead observing a larger region. HESS observed the Galactic Ridge (GR), defined 
as the region corresponding to longitude |£| < 0.8° and latitude |6| < 0.3°, with a cone of angle 
0.1° centered on the GC subtracted. 

This HESS observation leads to the dominant constraint on DM decays. Indeed HESS finds 
comparable total 7 fluxes from the GC and the GR regions, d^Qji/dE^ ~ Sd^cc/dE^, and we 
find comparable values of J: 

GR J for DM decays = {5.55, 20.5, 26.8} for the {isoT, NFW, Moore} profiles. (15) 

Therefore, since A^qr ~ 30Af2GC; observations of the GR set a bound on the DM life-time r 
which is about 10 times stronger than the GC bound. The blue dot-dashed lines in fig. 7 show 
our precise result. 

The GR HESS bound is close to the value of the DM life-time r suggested by PAMELA/ATIC. 
Such bound can presumably be improved by observing some region away from the Galactic 
Plane, to be chosen trying to avoid the astrophysical background. Unfortunately, astrophysi- 
cists are interested in astrophysical 'backgrounds', so that the GR region observed by HESS, 
that contains the galactic plane, is clearly not the optimal choice. 

4.3 HESS observation of Sagittarius Dwarf 

Dwarf spheroidals are among the most DM-dominated structures, so that they allow to search 
for 7 ray signals of DM annihilations with minimal astrophysical backgrounds. HESS observed 
Sagittarius Dwarf at distance d = 24kpc from us for a time Tobs = 11 h in a region with angular 
size Afl = 2 10^^, finding no 7 excess and setting the bound < 85 at about 3a. The 7 flux 
from DM decays is 



where Aes{E) ~ 10^ m^ is the effective area of HESS [16]. The integral is the total DM mass ^ 
contained in the volume observed by HESS, and is equal to ^ ~ 5 IO^Mq assuming a NFW 
or cored density profile [17] for Sagittarius Dwarf. Therefore, the bound on the DM life-time 
is: 




(16) 




(17) 



The resulting bound are plotted in fig.s 7 as blue dotted lines. 
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Figure 7: Assuming Dark Matter decays into e~^e~ (left column), ^'^ (middle column), r^r^ 
(right column) we compare the region favored by PAMELA (green band) and ATIC (red re- 
gion), as computed for min/max/med galactic propagation, with the bounds from HESS 
gamma observations of the Galactic Center (blue continuous curves), of the Galactic Ridge 
(dot-dashed curves) of Sagittarius Dwarf (blue dashed curves), and with radiowave observations 
of the Galactic Genter. We assumed NFW (upper row) and isothermal (lower row) DM density 
profiles. 
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4.4 Synchrotron radiation from the Galactic Center 



The produced by DM decays within the galactic magnetic field radiate synchrotron radiation. 
As the galactic magnetic fields are poorly known, in order to discuss which bounds are robust 
we recall the basic physics. One expects B ~ fiG that grows to higher values close to the 
galactic center. Such magnetic fields are intense enough that radiate an order unity fraction 
of their energy into synchrotron radiation, as other energy losses and diffusion are not the main 
phenomena. Therefore the total energy into synchrotron radiation can be robustly computed, 
and the value of the magnetic field only determines how it is distributed in the frequency 
spectrum. We recall that the synchrotron power radiated orthogonally to the magnetic 
field B by an with momentum p is: 

Thereby a larger magnetic field and heavier DM mean that synchrotron radiation extends up 
to higher energy. In order to put a robust bound, we consider the observation [18] of the GC 
in a region with angular size 4" at the lowest radio-wave frequency v = 0.408 GHz, that implies 
the bound S = {u dWsyn/diy)/{ATTrl) < 2 lO'^*^ erg/cm^ sec. This is the relevant bound, as the 
observed GC microwave spectrum is harder than what DM decays can produce. The region is 
large enough that we can neglect diffusion and advection (e^ falling in the central black hole). 
In this approximation the electron number density is given by ne{r,p) ~ E~y^ dE' Qe{r, E') 

with i^syn = e'^B'^p^ /^nm\ is the synchrotron energy loss of a in a turbulent magnetic field. 
So 

^^ = ^/ dVppNM (19) 



dv M 

where the integral extends over the observed volume and p is obtained from eq.s (18) as p 



A-Kmlv/B = 0.43 GeV(z//GHz)^/^(i?/mG)"^/^. Here, N^ip) is the number of electrons with 
energy above p produced in one DM decay, and can be often approximated as the total number 
of electrons produced in one DM decay. As anticipated at the beginning of this section, a lower 
B leads to a higher synchrotron fiux at the low frequency we consider. 

Observations are made at scales large enough that the dipolar structure of the GC can be 
neglected: we model the GC as in [9] in 'spherical-cow' approximation, obtaining the bounds 
plotted as red lines in fig.s 7. The radio-wave constraints turn out to be subdominant with 
respect to the 7-ray constraints previously computed. 



5 Theory 

In the DM annihilation case, model building is needed to invent models of DM that annihilate 
dominantly into leptons with cross section larger than the one suggested by cosmology. Key 
ingredients of such models are the Sommerfeld enhancement and possibly new vectors that 
decay into leptons [5, 6, 7]. 

DM decay allows to circumvent such model-building issues. It is very easy to invent models 
of DM that decays dominantly to leptons, and small couplings to leptons might be somehow 
connected to the small observed neutrino masses and baryon/lepton cosmological asymmetries. 
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One example is analogous to type-II see-saw [19]: DM might be a scalar SU(2)l triplet T 
with Lagrangian couplings M'^TT* + X^TLL + MXrpTHH, so that it dominantly decays into 
leptons if Xl ^ At- Alternatively, along the lines of type-I or type-III see-saw [19], one can 
introduce fermion singlets or triplets and an inert Higgs H' with couplings (for simplicity 
we consider the supersymmetric case and write the superpotential) MN'^/2 + XNLH' + eN^. 
The neutral component of H' could be DM and decay into LLL with a slow rate suppressed 
by (e/M)^. Within supersymmetry, one can build models of gravitino DM that decays into 
leptons, by assuming an appropriate sparticle spectrum and that LLE operators dominate 
i?-parity violation. DM might be some Froggatt-Nielsen-like field, and its slow decays into 
leptons might be related to the quasi-conserved U(l)£ flavor lepton numbers present in the 
lepton sector. 

We now focus on a more interesting idea. 

5.1 Non-annihilating DM and the Qdm/^b ratio 

A rather attractive scenario for decaying DM arises when one tries to relate the DM and the 
baryon energy densities in order to explain the ratio Qdu/^b ~ 5 [20]. 

We know that the amount of baryons in the Universe fi^ ~ 0.04 is determined solely by 
the cosmic baryon asymmetry UB/n^ ~ 6 x 10~^°. This is because the baryon-antibaryon 
annihilation cross section is so large, that virtually all antibaryons annihilate away, and only 
the contribution proportional to the asymmetry remains. This asymmetry can be dynami- 
cally generated after inflation. In contrast, we do not know if the DM density is determined 
by thermal freeze-out, by an asymmetry, or by something else. Thermal freeze-out needs a 
(Tf ~ 3 10~^^ cm^/ sec of electroweak size, suggesting a DM mass in the TeV range. If Qbu is 
determined by thermal freeze-out, its proximity to is just a fortuitous coincidence and is 
left unexplained. 

If instead ^dm ~ is not accidental, then the theoretical challenge is to define a consistent 
scenario in which the two energy densities are related. Since Qb is a result of an asymmetry, then 
relating the amount of DM to the amount of baryon matter can very well imply that Qbm is re- 
lated to the same asymmetry that determines Qb- Such a condition is straightforwardly realized 
if the asymmetry for the DM particles is fed in by the non-perturbative electroweak sphaleron 
transitions, that at temperatures much larger than the temperature T^, of the electroweak phase 
transition (EWPT) equilibrates the baryon, lepton and DM asymmetries. Implementing this 
condition implies the following requirements: 

1. DM must be (or must be a composite state of) a fermion, chiral (and thereby non-singlet) 
under the weak SU(2)i, and carrying an anomalous (quasi)-conserved quantum number 
B'. 

2. DM (or its constituents) must have an annihilation cross section much larger than elec- 
troweak (Tann ^ 3 10~^® cm^/ sec, to cusurc that ^dm is determined dominantly by the 
B' asymmetry. 

The first condition ensures that a global quantum number corresponding to a linear com- 
bination of B, L and B' has a weak anomaly, and thus DM carrying B' charge is produced in 
anomalous processes together with left-handed quarks and leptons [21, 22]. At temperatures 



14 



T ^ electroweak anomalous processes are in thermal equilibrium, and equilibrate the vari- 
ous asymmetries Fab = c^Ial = cb'Y^b' ~ O(10~^°). Here the Y^'s represent the difference 
in particle number densities n — n normalized to the entropy density s, e.g. Fab = {^b — ub)/ s. 
These are convenient quantities since they are conserved during the Universe thermal evolution. 

At T Mdm all particle masses can be neglected, and cl and cb' are order one coefficients, 
determined via chemical equilibrium conditions enforced by elementary reactions faster than 
the Universe expansion rate [23] . These coefficients can be computed in terms of the particle 
content, finding e.g. d = -28/51 in the SM and cl = -8/15 in the MSSM. 

At T -C Mdm, the B' asymmetry gets suppressed by a Boltzmann exponential factor 
q-Mj3m/t_ ]jgy feature of sphaleron transitions is that their rate gets suddenly suppressed 
at some temperature slightly below the critical temperature at which SU(2)i starts to be 
spontaneously broken. Thereby, if Mdm < T^, the B' asymmetry gets frozen at a value of 
0{Yab), while if instead Mdm > ^* it gets exponentially suppressed as Fa_b'/^ab ~ 

More in detail, the sphaleron processes relate the asymmetries of the various fermionic 
species with chiral electroweak interactions as follows. If B', B and L are the only quantum 
numbers involved then the relation is: 

where the order-one cl,b' coefficients are related to the cl,b' above in a simple way. The 
explicit numerical values of these coefficients depend also on the order of the finite temperature 
electroweak phase transition via the imposition or not of the weak isospin charge neutrality. 
In [24, 25] the dependence on the order of the electroweak phase transition was studied in two 
explicit models, and it was found that in all cases the coefficients remain of order one. The 
statistical function S is: 

-J I /o°° cosh~^ \\^^^ + -2^) for fermions , ^^^^ 

1 4^2 /o°° dx x^sinh"^ (^Vx^ + for bosons . 

with S'(O) = 1(2) for bosons (fermions) and S{z) ^ 12 [z /2TxY^'^e~^ at 2; ^ 1. We assumed the 
Standard Model fields to be relativistic and checked that this is a good approximation even for 
the top quark [25, 24]. The statistic function leads to the two limiting results: 

Y^B' ^ / 5(0) for Mdm <T, ...^ 

Fab 1 12 (MDM/27rT,)^/' e-^WT. for Mdm » T, ' ^ ^ 

Under the assumption that all antiparticles carrying B and B' charges are annihilated away we 
have Fab'/^ab = nB'/uB- The observed DM density 

f^DM _ Mdm nB' ^ ^ ^^g^ 
B Trip n B 

(where 1 GeV) can be reproduced for two possible values of the DM mass: 

i) Mdm ~ 5 GeV if Mdm ^ T^, times model dependent order one coefficients. 
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ii) Mdm ^ 8T^: ^ 2 TeV if Mdm ^ ^* , with a mild dependence on the model-dependent 
order unity coefficients. 

The ffist solution is well known [21] and not interesting for our purposes. The second solution 
matches the DM mass suggested by ATIC, in view of T^, ~ t; [22], where v ~ 250 GeV is the 
value of the electroweak breaking order parameter^. 

5.2 DM mass from strongly interacting dynamics 

The ideal DM candidate suggested by the PAMELA and ATIC anomalies, and compatible with 
direct DM searches, is a ~ 2 TeV particle that decays dominantly into leptons, and that has a 
neghgible couphng to the Z. 

If DM is an elementary particle, this scenario needs DM to be a chiral fermion with SU(2)i 
interactions, which is very problematic. Bounds from direct detection are violated. Further- 
more, a Yukawa coupling A of DM to the Higgs gives the desired DM mass Mdm ~ Af ~ 2 TeV 
if A ~ 47r is non-perturbative, hinting to a dynamically generated mass associated to some new 
strongly interacting dynamics [27, 22, 24, 25]. This assumption also solves the problem with 
direct detection bounds, which are satisfied if DM is a composite SU(2)2,-singlet state, made of 
elementary fermions charged under SU(2)i,. 

This can be realized by introducing a strongly-interacting ad-hoc 'hidden' gauge group. A 
more interesting identification, namely technicolor, is suggested by the the proximity of the DM 
mass indicated by ATIC, Mdm ~ 2 TeV, with the Anv scale at which strong dynamics might 
naturally generate the breaking of the electroweak symmetry. In such a scenario, DM would 
be the lightest (quasi)-stable composite state carrying a B' charge of a theory of dynamical 
electroweak breaking featuring a spectrum of technibaryons {B') and technipions (11). 

Let us elaborate more quantitatively on this numerical connection. The DM mass can be 
approximated as m^' = Mdm ~ ^qAtc where uq is the number of techniquarks Q bounded 
into B' and Atc is the constituent mass, so that MDu/mp ~ tt-qAtc/SAqcd- Denoting by 



f-K (Fyi) the (techni)pion decay constant, we have F^/ = y D/S Atc/Aqcd where Dq is the 
dimension of the constituent fermions representation {D = 3 in QCD)^. Finally, the electroweak 
breaking order parameter is obtained as v"^ = Nj:,F^, from the sum of the contribution of the 
Nd electroweak techni-doublets. Putting all together yields the estimate: 



where the numerical value corresponds to the smallest number of constituents and of techni- 
quarks hq = Dq = 2 and A''^ = 1. We conclude that naive rescaling of QCD yields a value of 
Mdm right in the ballpark suggested by ATIC.^ 

^More precisely, for a Higgs mass mh = 120 (300) GeV, Ref. [26] estimates « 130 (200) GeV, where the 
larger values arise because of the larger values Higgs self coupling. For the large masses that are typical 
of composite Higgs models, the self coupling is in principle calculable and generally large [25], so that taking 

^ w is not unreasonable. 

^The large- counting relevant for a generic extension of technicolor type can be found in Appendix F of [28] 
together with a general introduction to recent models of dynamical electroweak symmetry breaking. 

^ Besides to the possibility indicated above, a dynamical origin of the breaking of the electroweak symmetry 
can lead to several other interesting DM candidates (see [28] for a list of relevant references). 





V 



m„ = 2.2 TeV 



(24) 
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5.3 Phenomenological constraints on techni-DM 

It is worth mentioning that models of dynamical breaking of the electroweak symmetry do 
support the possibility of generating the experimentally observed baryon (and possibly also the 
technibaryon/DM) asymmetry of the Universe directly at the electroweak phase transition [29]. 
Electroweak baryogenesis [30] is, however, impossible in the Standard Model [31] (see [32] for 
a review on this topic). 

Bounds from direct detection experiments such as CDMS and Xenon suggest that a single 
component DM candidate can have, at most, a tiny coupling with the photon [33], the Z [34, 
35] and the gluons. We therefore assume that our technibaryon B' is a singlet of 811(2)/, 
with zero hypercharge (and thus also electrically neutral). Being composed of fermions with 
electroweak interactions, analogously to the neutron, the DM particles will have electro (weak) 
magnetic dipoles and other similar form factors, that can be parameterized at low energy by 
non-renormalizable effective operators suppressed by powers of the technicolor scale Atc- For 
example, if the B' particle is a composite scalar, the relevant operators are suppressed at least by 
two powers of Atc [36], and it has been checked in [24] that in this case the expected interaction 
rate remains below the present bounds [37]. An interesting collection of other techno-cosmology 
estimates can be found in [27]. 



5.4 DM lifetime and decay modes 

According to [38] the sphaleron contribution to the techni-baryon decay rate is negligible be- 
cause exponentially suppressed, unless the techni-baryon is heavier than several TeV. 

Grand unified theories (GUTs) suggest that the baryon number B is violated by dimension-6 
operators suppressed by the GUT scale Mqut ~ 2 ■ 10^^ GeV, yielding a proton life-time [39] 

^(p ^ ;r°e+) ~ ~ 10^1 sec. (25) 

If B' is similarly violated at the same high scale Mqut, our DM techni-baryon would decay 
with life-time 

--^-10"-e, (26) 

which falls in the ball-park required by the phenomenological analysis above. Models of unifi- 
cation of the Standard Model couplings in the presence of a dynamical electroweak symmetry 
breaking mechanism have been recently explored [40, 41]. Interestingly, the scale of unification 
suggested by the phenomenological analysis emerges quite naturally [41]. 

Low energy DM and nucleon (quasi)-stability imply that, in the primeval Universe, at 
temperatures T < Mgut perturbative violation of the B' and B global charges is strongly 
suppressed. Since this temperature is presumably larger than the reheating temperature, it 
is unlikely that Qb and I^dm result directly from an asymmetry generated in B' or B. More 
likely, the initial seed yielding Qbm and Qb could be an initial asymmetry in lepton number 
L that, much along the lines of well studied leptogenesis scenarios [42], feeds the B and B' 
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asymmetries through the sphaleron effects.^ Indeed, it has been shown that it is possible to 
embed seesaw-types of scenarios in theories of dynamical symmetry breaking, while keeping the 
scale of the L- violating Majorana masses as low as ~ lO^TeV [46]. 

Assuming that techni-baryon DM decay is dominantly due to effective four-fermion oper- 
ators, its decay modes significantly depend on the technicolor gauge group. In the following 
F generically denotes any SM fermion, quark or lepton, possibly allowed by the Lorentz and 
gauge symmetries of the theory. 

• If the technicolor group is SU(3), the situation is analogous to ordinary QCD: DM is a 
fermionic QQQ state, and effective QQQF operators gives DM — > decays. This 
leads to hard leptons, but together with an excess of p, from the n~ — sc decay (in view 
of n~ ~ W£). This is therefore incompatible with PAMELA p data, unless astrophysical 
propagation uncertainties are relaxed. 

• If the technicolor group is SU(4) the situation is typically worse: DM is a bosonic QQQQ 
state, and effective QQQQ operators lead to its decay into techni-pions. 

• Finally, if the technicolor gauge group is SU(2), DM is a bosonic QQ state, (as put 
forward in [24]), and effective QQFF operators lead to DM decays into two F. Since 
the fundamental representation of SU(2) is pseudoreal, one actually gets an interesting 
dynamics analyzed in detail in [24]. Here the technibaryon is a pseudo-Goldstone boson 
of the underlying gauge theory. 

An SU(2) technicolor model compatible with the desired features is obtained assuming that the 
left component of the Dirac field Q has zero hypercharge and is a doublet under SU{2)l, so that 
DM is a scalar QQ with no weak interactions, and the four-fermion operator {QQ)d ^^{F'-) ^F) 
allows it to decay. Such operator is possible for both SM leptons and quarks, so that the DM 
branching ratios into and qq is a free parameter. 

The use of higher dimensional representations for the techniquarks transforming under a 
given technicolor gauge group has opened new possibilities [28]. Interestingly, the technicolor 
model with the lowest possible value for the naive S parameter (measuring deviations from 
electroweak precision data) can be constructed by combining exactly two Dirac flavors in the 
fundamental of SU(2) with one Dirac in the adjoint representation (uncharged with respect to 
the SM) [24]. The model has also the lowest number of fermions and is compatible with near 
conformal dynamics, which is an important ingredient for reducing the tension with constraints 
stemming from unobserved flavor changing neutral currents. 

6 Conclusions 

We explored the interpretation of the excess of cosmic rays observed by PAMELA and ATIC 
in terms of decays of Dark Matter with mass Mdm- 

'*In Minimal Walking Technicolor [43, 44], one additional (techni-singlet) SU( 2) -doublet must be introduced 
to cancel the odd-number-of-doublets anomaly [45]. An asymmetry in the L' global charge associated with 
these new states can also serve as a seed for the B and B' asymmetries. 
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On the phenomenological side, we found that any DM decay modes involving hard leptons 
can fit the PAMELA excess for M^m ^ 200 GeV, while decay modes that give soft leptons with 
energy E <^ Mdm together with p need M^m ^ 10 TeV. The peak present around 800 GeV 
in the e"*" + e~ ATIC spectrum can be reproduced if Mdm ~ 2 TeV and for DM decays into 
H'^fi~ , T^T~ (characteristic of boson DM) and into W^t^ where £ is any lepton (characteristic 
of fermion DM). However, the latter decays are disfavored by the PAMELA non-observation of 
an excess in the p CR spectrum. Fermionic DM is allowed if it decays into i~^i~i>. 

DM annihilations compatible with the PAMELA excess tend to give a fiux of 7 rays above 
the HESS [17, 47, 48] bounds and a fiux of synchrotron radiation above the Davies [18] bounds, 
unless the DM density profile p{r) is significantly less steep as r — >■ than the NEW profile. 
Whether or not this should be considered a problem depends on the the reliability of A^-body 
simulations, and of their extrapolation to scales below 100 pc. 

Anyhow, this issue is not present for DM decays, because they give DM signals proportional 
to p rather than to p^. The HESS observation of the Galactic Ridge region gives the dominant 
constraint, which is just below the needed sensitivity, and can be improved by performing an 
observation optimized for DM rather than for astrophysics. 

On the theoretical side, we found that it is easy to invent decaying DM models compat- 
ible with the above features suggested by data, just because model-building is less severely 
constrained than in the case of DM annihilations. 

More interestingly, if the DM density is due to an asymmetry analogous to baryon number 
broken at low energy by weak anomalies, then the value of the DM mass suggested by the 
e"*" + e~ peak, Mdm ~ 2 TeV, is also the one that naturally gives the observed ratio between 
the dark and 6aryonic matter densities VI-dm/^b ~ 5. 

Indeed, in such a case electroweak sphalerons keep tzdm ~ down to T>Mdm; below 
this temperature the DM density gets Boltzmann-suppressed, until sphalerons freeze-out at the 
temperature ~ 200 GeV. As a result, the correct DM density is obtained for M ~ 8T=k, which 
coincides with the mass suggested by the ATIC peak. Technicolor extensions of the SM suggest 
a similar value for the mass of the lightest (quasi)-stable state. Technicolor is a new strongly 
coupled gauge dynamics which breaks the electroweak symmetry dynamically and solves the 
hierarchy problem without an elementary Higgs. These models feature technibaryons which are 
quasi-stable technihadrons analogous to the proton in QCD, that possess, quite surprisingly, 
many of the required properties advocated by our phenomenological analysis. 

While these phenomenological coincidences are indeed suggestive, from the theoretical point 
of view, bound states of some new strong dynamics decaying predominantly into leptons require 
some specific construction, that seems to privilege an SU(2) technicolor gauge theory with a 
technimatter structure envisioned in [24]. 
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